Direct measurement of shear-induced cross-correlations of Brownian motion 
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Shear-induced cross-correlations of particle fluctuations perpendicular and along stream-lines are 
investigated experimentally and theoretically. Direct measurements of the Brownian motion of 
micron-sized beads, held by optical tweezers in a shear-flow cell, show a strong time-asymmetry 
in the cross-correlation, which is caused by the non-normal amplification of fluctuations. Comple- 
mentary measurements on the single particle probability distribution substantiate this behavior and 
both results are consistent with a Langevin model. In addition, a shear-induced anti-correlation be- 
tween orthogonal random-displacements of two trapped and hydrodynamically interacting particles 
is detected, having one or two extrema in time, depending on the positions of the particles. 

PACS numbers: 5.40.Jc, 82.70.-y, 47.15. G-, 87.80.Cc 



The Brownian motion of particles in fluids and their 
hydrodynamic interactions are of central importance in 
chemical and biologicalphysics as well as in material sci- 
ence and engineering [l|, 0, 0, [3- However our under- 
standing of the dynamics of particles in flows is still far 
from complete. Direct observations of particles at the 
mesoscale substantially contribute to our understanding 
of their dynamics. At this scale optical tweezers are a 
powerful experimental technique [5| with a number of 
innovative applications. They include the detection of 
anti-correlations between hydrodynamically interacting 
Brownian particles , propagation of hydrodynamic in- 
teractions [7| , short-time inertial response of viscoelastic 
fluids [1], two-point microrheology @, anomalous vibra- 
tional dispersion [l^ and particle sorting techniques [ll| . 

Neutral colloidal particles moving relatively to each 
other interact via the fluid and these hydrodynamic in- 
teractions decay with the particle distance [3. In shear 
flow little is known about the dynamics of Brownian par- 
ticle motion and the hydrodynamic interaction effects in 
spite of their fundamental relevance and importance in 
applications in microfluidics, Taylor dispersion and 
in fluid mixing [1, [l^ . In time dependent fields and in 
shear flow surprising deterministic particle dynamics may 
be induced by hydrodynamic interactions [l4| . For poly- 
mers it is the interplay of shear flow and fluctuations 
which leads already at low Reynolds numbers, to rich 
dynamics , the so-called molecular individualism , 
causing elastic turbulence even in diluted polymer solu- 
tions [17[ and spectacular mixing behavior [l3|. 

It is the contribution (u • V)u to the Navier-Stokes 
equation which causes interesting transie nt p henomena 
in shear flows near the onset of turbulence [1^ , as well as 
amplifications of fluctuations and their cross-correlations 
along and perpendicular to straight streamlines [l^, [2^ . 
A cross-correlation is also expected between orthogonal 
particle-fluctuations in the shear plane, because random 
jumps of a particle between neighboring streamlines of 
different velocity lead to a change of the particle's ve- 
locity and displacement along the streamlines, similar 
as via fluctuations. In some parameter ranges, iner- 



tia effects may become important (2l|, |22|, |23[. Cross- 
correlations between perpendicular fluid-velocity fluctu- 
ations and perpendicular fluctuations of particles are ex- 
pected to be strongly asymmetric in time [l^, [2^, [23 |. 
In dynamic light-scattering experiments certain aspects 
of these shear-induced cross-correlations were observed 
indirectly [25| . but a direct measurement and character- 
ization of related particle fluctuations is missing. 

Here we investigate in a linear shear flow the fluctua- 
tions of a single particle in a potential minimum and of 
two hydrodynamically interacting particles trapped by 
two neighboring potentials. We use a special shear-flow 
cell, where one or two micron sized beads are held at its 
center by optical tweezers. The time-asymmetry of shear- 
induced cross-correlations were determined directly by 
measuring the particle's positional fluctuations. In ad- 
dition the probability distribution of a single particle in 
a trap was measured, which can be also calculated in 
terms of a Langevin model, similar to the correlations. 
Both the probability distribution and the correlation can 
be fitted by using the same value of the shear rate, which 
altogether gives a consistent picture of not yet directly 
observed shear-induced cross-correlations of particle fiuc- 
tuations. 

By a dual beam optical tweezer-setup, composed of two 
solid state lasers and an oil immersion objective with a 
numerical aperture of 1.4, two harmonic potentials are 
generated in an inverted microscope (Nikon TE 2000-S) 
to capture uncharged polystyrene beads (Duke Scientific 
Corporation, R0300) with a diameter of 3/xm in a fiow 
of distilled water. The beads were observed with a high 
speed camera (IDT, X-Stream, XS-5) of 15 kHz and their 
positions were determined with a correlation tracking al- 
gorithm with a spatial resolution of ±4 nm [2^ . To avoid 
any interference of the two potentials at small distances 
and to maximize the hydrodynamic interaction effect be- 
tween the two particles, they where held at a distance of 
d = 4/iTO. In a microfiuidic device with two counter fiows, 
as shown in Fig. (TJ a linear shear gradient with a van- 
ishing mean velocity was generated at the center of the 
cell, as experimentally verified by micro-PIV (see inset in 
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Fig. [T]) . The design of the flow chamber was optimized 
by numerical simulations of the incompressible Navier- 
Stokes equation (Multiphysics 3.4, Comsol AB, Stock- 
holm, Sweden). The small width in x direction of the 
center piece of the chamber was chosen in order to mini- 
mize flow in y direction. The curved form of the bound- 
aries was found to suppress vortices. The channel was 
manufactured by standard soft lithographic techniques 
and the flow was driven by gravitational potential differ- 
ence. The distance of the beads from the wall was always 
larger than lO^m in the z direction and 25fj,m in the xy 
plane. Hence boundary effects on the bead fluctuations 
could be excluded within our experimental resolution, as 
verified by measurements without flow, which were in 
very good agreement with previous results (see e.g. Ref. 
0). However, with flow the experimental noise becomes 
larger, especially at longer correlation times. The value 
of the shear rate 7 has been extracted from the fits of the 
correlation data and by particle tracking methods. The 
highest shear rate that did not lead to an escape of the 
particles from the traps was 7 ~ 50s~^ [27| . 
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FIG. 1: A cell cross-section (lithographic mask) is shown with 
150/im depth in the z direction having opposite flow directions 
in its upper and lower channel and a linear shear profile at 
its center: see inset with PIV data (A) and linear fit. One 
or two particles at a distance d = Afim were held by optical 
tweezers in the center of the linear velocity profile u(y). 



One or two Brownian particles with coordinates = 
(xi, j/i, Zi) [i = 1, 2) are held in a linear shear flow \i{yi) — 
jyiSx by forces = k{pi — r^) close to the minima Pi of 
two harmonic potentials Vi = |(pi— r,)^ (spring constant 
k). The over-damped particle motion is described by a 
Langcvin equation 



r, = u(r,) + (f/ + f/) 



(1) 



The mobility matrix accounts for the Stokes friction 
and the hydrodynamic interactions between them. Here 



we use the Oseen approximation 
Hii = H22 ^ ' 



H12 = H 
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with the Stokes friction coefflcient C = Gmja of a point 
particle of effective hydrodynamic radius a in a fluid of 
viscosity 77 and the unity matrix E. ri2 = ri — r2 is 
the bead distance and ri2 is its norm, t = C^/k the par- 
ticle relaxation time in the potential and W = jr the 
Weissenberg number. The Brownian particle motion is 
driven by the stochastic forces if {t) in Eq. ([T]), for which 
we assume vanishing mean values and correlation times: 
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At first we investigate the Brownian motion of a sin- 
gle trapped particle in shear flow. Its autocorrelation 
along the flow direction, (a;(0)a;(0)) = ^(1 + 
depends on W, but along the perpendicular direction, 
(y(t)y(0)) = exp(— |:), it docs not depend on 

W [2^. In a quiescent fluid cross-correlations between 
particle displacements in orthogonal directions vanish: 
{x{t)y{t')) = 0. But shear flow causes in the shear plane 
finite cross-correlations i2L. |23i . i28i] , which are asymmet- 
ric with respect to t ^ —t [24| : 
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The algebraic prefactor in Eq. ([6]) illustrates that a fiuc- 
tuation y{0) 7^ of a particle is carried away by the 
fiow in the X direction before the initial displacement 
y(0) relaxes. This leads, during an initial period shorter 
than the relaxation time r, to a growth of {x(t)y{0)), 
while the expression in Eq. ^ decays monotonically. As 
shown in Fig. [21 the predicted elementary signatures for 
the shear-induced cross-correlations, cf. Eq. ([6]), are in 
agreement with our experimental data (triangles). Here 
{x{t)y{Q)) takes its maximum roughly at t w 0.009s, cor- 
responding via Eq. ([S]) to a particle's relaxation time 
T w 0.018s. Also the initial decay of {x{0)y{t)) (squares 
in Fig. [2]) agrees with our model, cf. Eq. ([7]). The addi- 
tionally observed minimum is possibly caused by a slight 
inclination of the laser beam or it is a reminiscent of a 
long wave length oscillation due to the limited number 
of samples taken For fluid- velocity fluctuations in 

orthogonal directions in the shear plane a similar sig- 
nature as in Ec^. ^ has been found [l^. According to 
Eq. ([6] and Eq.[7| one obtains the normalized ratios of the 
static cross-correlations: (.■r(0)y(0))/(y(0)j/(0)) — W/2 
and {x{0)y{0))/{x{0)x{0)) = M- From the flts, 

as indicated by the red and blue line in Fig. [H we ob- 
tain (x(0)j/(0))/(a;(0).T(0)) 0.26, which corresponds to 
a Weissenberg number W w 0.62. 
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FIG. 2: Shear-induced cross-correlations between orthogonal 
random displacements: {x{t)y{0)} (triangles) and {x{0)y{t)) 
(squares). Lines are fits according to Eq. ((6| and Eq. ((7)l and 
open circles represent {x{0)y{t)) in the absence of flow. 



The probability distribution of a Brownian particle in 
a harmonic potential and exposed to a linear shear flow 
has an elliptical cross-section as shown by the particle's 
position in Fig. [3] but it has circular symmetry in the ab- 
sence of flow. The angle enclosed by the major axis of 



O.lr 






0.1 




E 
>, 















Y =0 


-0.1 



b) 




Y = 49s" 



-0.1 



FIG. 3: The particle's distribution in the shear plane is shown 
in a) without flow and for a shear flow with 7 = 49s ~^ in b). 
The angle between the major and the x axis is </> « 38° and 
the ratio between the two principal axes is i? ~ 0.75. 

the particle's probability distribution and the x axis, as 
well as the ratio R between the lengths of the two princi- 
pal axes, in the shear plane, depend on the Wcisscnberg 
number W as follows 1241: 



tan 6 ^ — 
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(8) 
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Using W 0.62 as determined above, one obtains via 
Eq. dS]) the angle 4> « 37° and via Eq. (O the ratio R « 
0.72. Within errors this is consistent with the angle (p « 
38° and the ratio R w 0.75 obtained from the measured 
particle's distribution shown in Fig. [31 



For two particles, each trapped in a potential minimum 
in shear flow, we investigated the correlations between 
their random displacements for two different configura- 
tions: with the connection vector pi2 = Pi — P2 parallel 
to the flow direction as in Fig. [T]or perpendicular to it. 

For Brownian displacements of the two distinct parti- 
cles along the same direction the quantities {xi{t), Xj{0)) 
and {yi{t)_jjj{Q)) describe anti-correlations for i ^ j (see 
e.g. Ref. 0). The shear-induced corrections for both are 
of the order of as described in more detail in Ref. [U . 
For random displacements of distinct particles, but along 
orthogonal directions, one only finds correlations in the 
presence of shear flow. With the abbreviations 



Ai,3 = l±2/.(, A2,4 = 1±M7 
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and the connection vector pi2 parallel to the flow two of 
the anti cross-correlations in the shear plane are [l^l 

(xi(0y2(0)) = ^^re-^^^/^ + e-^^*/^ 
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The cross-correlation {xi{t)y2{0)) (triangles) in Fig. [4] 
and the fit (blue line) show a pronounced minimum at 
about the particles relaxation time t « r. 
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FIG. 4: Correlations {xi(t)y2{0)) (triangles) and (xi (0)1/2 (t)) 
(squares) between random displacements of two particles. 
Colored lines are fits according to Eq. and Eq. (|12p . Cir- 
cles represent the same correlations in the absence of flow. 

With a connection vector pi2 perpendicular to the flow 
lines we obtain a cross-correlation {x2(t)yi (0) in the limit 
of small values of ^, 



t 
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2k 2 

which exhibits in contrast to Eq. pip two cxtrema 
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Shear-induced cross-correlations between random dis- 
placements of a single particle in a potential were calcu- 
lated and measured here for the first time, to the best 
of our knowledge, cf. Fig. [51 At approximately half of 
the particle's relaxation time t the correlation function 
in Eq. (jS]) exhibits with its maximum a typical signature 
of Brownian motion in shear flow, caused by the rota- 
tional part of the shear flow as well as the non-normal 
property of the linearized Navier-Stokes equation. Si- 
multaneously, for a particle in a harmonic potential and 
shear flow an elliptical probability distribution was mea- 
sured. Both independent measurements are described by 
a Langevin model for the same value of the Weissenberg 
number, which confirms the validity of our approach to 
shear-flow effects on the Brownian particle dynamics. 

Theoretically, shear-induced correlations between per- 
pendicular fluid-velocity fluctuations have been investi- 
gated before [l^, [l^l • Those are traced back to the non- 
normal property of the linearized Navier-Stokes equation 
and they are important for the stability of shear flow 
and the onset of turbulence. The cross-correlations be- 
tween these velocity fluctuations are based on the same 
mechanism as discussed here and they exhibit similar ex- 
trema as our experimental and analytical results. 

Stochastic forces on a suspended particle are caused 



by velocity fluctuations of the surrounding fluid. Usu- 
ally, they are assumed to be isotropic in related Langevin 
models with uncorrelated perpendicular components. 
However, cross-correlations of the velocity fluctuations in 
shear flow, as discussed in Refs. will modify the 

cross-correlations between orthogonal particle displace- 
ments, as investigated here, but the related additional 
contributions to the particle displacement correlations 
are expected to be considerably smaller than the effects 
of isotropic random forces [IJ]. It is, however, an inter- 
esting and challenging future issue to separate these two 
non-equilibrium effects in experiments. 

For two hydrodynamically interacting particles, each 
captured by an optical tweezer at the center of the shear 
flow, we find shear-induced anti-correlations between or- 
thogonal particle displacements with one cxtrcmum if the 
vector connecting the mean particle positions is parallel 
to the streamlines and two extrema, if the connection 
vector is perpendicular to the flow lines. These prop- 
erties may be relevant for further understanding of the 
dynamics of polymer models in shear flow. 

This work was supported by the German science foun- 
dation via the priority program on micro- and nanoflu- 
idics SPP 1164 and the graduate school GRK 1276. 
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